The fate of pion condensation in quark matter: 
from the chiral limit to the physical pion mass 
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We study aspects of the pion condensation in two-flavor neutral quark matter using the Nambu- 
Jona-Lasinio model of QCD at finite density. We investigate the role of electric charge neutrality, 
and explicit symmetry breaking via quark mass, both of which control the onset of the charged 
pion (n c ) condensation. We show that the equality between the electric chemical potential and the 
in-medium pion mass, fi e = M w - , as a threshold, persists even for a composite pion system in the 
medium, provided the transition to the pion condensed phase is of the second order. Moreover, we 
find that the pion condensate in neutral quark matter is extremely fragile to the symmetry breaking 
effect via a current quark mass m, and is ruled out for m larger than the order of 10 keV. 
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I. INTRODUCTION 



In the early 1970s Migdal suggested thepossibility of 
pion condensation in a nuclear medium [l|, Since then, 
many efforts have been made in order to clarify the in- 
mcdium pion properties affected by the pion-nucleon in- 
teraction Q, because a better insight of such conden- 
sation phenomena would yield important advances not 
only in subnuclear physics such as that of pionic atoms 
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4 1 , but also in the physics of neutron stars , supernovas 
6 j, and the heavy ion collisions Q. From the results of 
the experiments concerning the repulsive ttN interaction, 
the simplest possibility, the S-w&ve pion condensation is 
highly unlikely to be realized in nature. In contrast, the 
possibility of P-wave pion condensation has been widely 
argued [HfilH , even though it still remains as a matter of 
debate if this occurs at several times the ordinary nuclear 
density. Most of the studies about this issue performed so 
far are mainly concerned with the role of the coupling be- 
tween pions and baryons in the nuclear medium Q , con- 
sidering the pion itself as an elementary object. However, 
since the pion is considered to be the Nambu-Goldstonc 
boson associated with chiral symmetry breaking, its in- 
ternal structure and mass may also be sensitive to the 
modification of the QCD vacuum itself in the finite den- 
sity environment [To| ; the finite baryon density, and even 
the isospin density arising from the neutrality condition, 
can modify the structure of the QCD ground state, and it 
can in turn produce significant modifications of the pion 
properties in the medium. 
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In this article we revisit the possibility of S-w&ve 
charged pion (7r c ) condensation starting from a micro- 
scopic model which is built with quarks as the con- 
stituents of pions, and which exhibits chiral symmetry 
restoration at the finite quark chemical potential [i or 
temperature T. This strategy also enables us to present 
a unified view on the possible crossover from a Bosc- 
Einstcin condensate (BEC) of ir~ condensate to a BCS- 
typc (uj^d) su pcr fluid (BCS-BEC crossover) at finite 
isospin density This possibility was first speculated 
by applying the chiral Lagrangian at low \ij and the 
perturbative QCD result at high \ij [l2j|; since then, 
the QCD phase structure at finite isospin density has 
been widely investigated using the chiral effective models 
[SOiEl , a lattice-based model [ll| , the random matrix 
model [l7| . and the lattice QCD simulations [lei [l9l. |20T] . 

Since \i acts as the external field producing a mismatch 
in u-d Fermi surfaces [HI, [2l| , our problem also closely 
intersects with the problem of the BCS-BEC crossover 
with a density imbalance in the cold atomic systems [22| . 
Finally, it should be noted that in Rcfs. [23|, [24| some 
window was found for the S'-wave n c condensate on the 
finite /z axis in the Nambu-Jona-Lasinio (NJL) model. 
This was called the "gapless pion condensate" since the 
constituent quark presents a gapless dispersion relation, 
due to u-d Fermi surface mismatch. Our purpose here 
includes clarification of how this can be understood con- 
sistently with the number of negative results on the S- 
wave pion condensate obtained in the past. To this end, 
we derive the appropriate criterion for 7r c condensation, 
and based on this criterion, we show that the 7r c conden- 
sation found in [23j is extremely fragile to explicit chiral 
symmetry breaking due to a finite current quark mass. 

The purpose of this article is twofold. First, we inves- 
tigate the conditions for the onset of 7r c condensation 
at finite density using the NJL model of QCD; in this 
model pions are not elementary fields, but they are de- 
scribed as the bound states of quarks. We find that even 
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in such a composite pion model with dynamical quarks, 
the threshold for 7r c condensation for noninteracting ele- 
mentary pion gas, /i e = M w - (— [i e = M^+) for positive 
(negative) fi e (with fj, e being the electric chemical poten- 
tial), holds as well. This enables us to elucidate the effect 
of the current quark mass in the neutral ground state, 
which will be summarized by drawing the phase diagram 
in the (/i, m w ) plane, where m^ denotes the vacuum pion 
mass. Furthermore wc investigate the effect of fi e at the 
physical value of the current quark mass drawing a phase 
diagram in the (/z, /Lt e ) plane. Based on these analyses, we 
conclude that 7r c condensation is forbidden in realistic 
neutral quark matter. 



II. THE MODEL 

We study two-flavor quark matter at a finite chemical 
potential within the NJL model. The Lagrangian of the 
model is given by [25[ 

£ = eii'y^d' 1 + pt e 7o)e + $ (ij^ + p,j - m) ip 



+G 



(tpip) 2 + (ipi-ysrip) 



(1) 



Here e denotes the electron field, and ip is the quark 
spinor with Dirac, color and flavor indices (implicitly 
summed), to = m u = ma is the bare quark mass and G is 
the coupling constant. fi e is the electric charge chemical 
potential needed to keep the system electrically neutral 
|25j | . while /z/ serves as the isospin chemical potential in 
the hadron sector, /!/ = — \i e since Q = + I3. The 
quark chemical potential matrix \x is defined in flavor- 
color space as /t = diag(/i — §/i e i£* + gMe) ® l c , where 
l c denotes the identity matrix in color space, /i is the 
quark chemical potential related to the conserved baryon 
number. By performing the mean field approximation, we 
examine the possibility that the ground state develops 
condensation in the a = G(ipip) and/or tv = G(ipi"/5Ti/j) 
channels, where t = {tx,T2,ts\ denotes the Pauli ma- 
trices. 1 We find that (773) is always zero, as there is 
no driving force, so wc omit it in the following. We 
use the notation, 



M 



2cr and N = 20rf 



In the numerical analyses, we fix A = 651 MeV and 
G = 2.12/A 2 so that the model reproduces / w = 92 MeV, 
(uu) = -(250 MeV) 3 , and = 139 MeV in the vacuum 
with m = 5.5 MeV. 



III. THE SYMMETRY STRUCTURE 

The global symmetry of the model at m = /i e = is 
the SUl(2) x SL/r(2) chiral symmetry. This symmetry 



1 In this study we restrict ourselves to the homogeneous ground 
state, although some interesting possibilities of spatially modu- 
lated condensate are argued in Ref. |26H . 



will be broken spontaneously to ST/ (2) by the emergence 
of a nonzero condensate v 'o~ 2 + ir 2 . When the finite quark 
mass to is turned on, the original chiral symmetry is ex- 
plicitly broken to the diagonal subgroup SUi,+r(2). This 
forces the chiral condensate to the a direction (vacuum 
alignment), and the pions as Goldstone bosons acquire a 
mass, but they remain degenerate in the isospin triplet. 
The effect of a small /ij is then to split the masses of the 
three pions, since the diagonal ST/l+r^) is broken ex- 
plicitly to U T3 (1), which denotes the rotation about the 
isospin third axis. When 7Ti and/or 112 condensates (or 
simply N) become nonzero, this U T3 (1) symmetry gets 
spontaneously broken and we expect one massless Gold- 
stone boson to appear. It is interesting to note that when 
the discrete parity is considered together, there remains 
one discrete symmetry even in the pion condensed phase. 
The parity operation acts as it — > Pit = —ir. On the other 
hand, Z 2 C U T3 (l) defined by U G Z 2 & U = e'7 T3 , also 
flips the sign of the 7r c condensate. Hence, the combined 
transformation, a rotated parity P' = P ■ Z2, remains un- 
broken. The new symmetry breaking pattern of the 7r c 

condensate can be specified as U T3 (1) x P ^4° P'. The 
new parity P' is analogous to what is discussed in [27l |. 
Finally, [i breaks the charge conjugation, causing stress 
which produces a mismatch in u and d densities that are 
equal at small (large) /i {pi). 



IV. PION CONDENSATION AT FINITE pj AND 
AT FINITE (i 

Before discussing tt c condensation in the medium, it 
is instructive to review the classical arguments at fx = 
[T3| ■ The shift of the vacuum energy caused by the pion 
fields at the lowest nonvanishing order in the fields tt is 
60, = (ml + /i 2 )7r^/4 + (ml - hI)tt + tt-/2 . From this we 
can infer that in the chiral limit (m w = 0) an infinitesi- 
mally small value of fij favors condensation, (7r + 7r_) 7^ 0. 
In the general case of m T 7^ and fj, e 7^ 0, the condition 
|/i e | > m-x needs to be fulfilled so that the 7r c condensa- 
tion turns out to be energetically favored. In this kind 
of approach, using chiral perturbation theory, the chi- 
ral SU (2) multiplets are elementary fields whose internal 
structures do not suffer from any quantum corrections. 
On the other hand, the same threshold at fx = is ob- 
tained also in the NJL- type model [Til ] . 

We now consider neutral quark matter at 11 7^ and 
T = 0. We are interested in the relation between the 
threshold of tt c condensation at finite density and the 
in-medium pion masses in the neutral ground state. We 
have shown in Ref. [25| that at the physical point m — 5.5 
MeV, there is no room for tt c condensation in the neutral 
phase. Similar conclusion is also obtained in [24j . This 
picture changes if we lower the current quark mass. We 
discuss this point in detail later. For our present pur- 
pose it is enough to state that we need a current quark 
mass of the order of 10 kcV. In Fig. Q] we plot M and 
V in the neutral phase as a function of fi for to = 10 
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keV. In this figure M n o, M v ± denote the in-medium pion 
masses defined by the poles of the pion propagators in the 
rest frame, computed in the random phase approximation 
(RPA) to the Bethe-Salpeter (BS) equation. For charged 
pion sector in the phase of a vanishing pion condensate 
N = 0, the pole condition is r „.+„.- (Qq; M, Me, T) = 
with r 7r+7r -(<5o;A 1 ^e,7 1 ) = - N c QlJ(Q ; fi u , Md,T) 
being the inverse ir + propagator, Qq = lu — /i e , N c be- 
ing the number of colors, and J being the polarization 
defined by 



j(<3 ;m 

V ( 9° 



-21 E 



2Qo 



f F (E-tfi u ) + 
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2Qo 



E 2 - 

f F (E- 



tUd) 



1 



(2) 

where Jf(x) is the Fermi distribution. We note that the 
polarization function is not symmetric under Qq — > — Qq 
if there is a finite isospin density in the system, (ip'Tsip) ^ 
0; this asymmetry is responsible for the mass splitting of 
M w + and M„- . The excitation gaps for ir + and ir~ are 
accordingly given by (M^+ + Me) and (M T - — Me): cor- 
responding to positive and negative solutions to the BS 
equation in w. From Fig. [T] we notice that the transition 
to the pion condensed phase is of second order and it 
occurs at the point where M T - = Me- 

The numerical results shown in Fig. [T] can be under- 
stood by a closer inspection to the BS equation and the 
gap equation for pion condensate. Although the argu- 
ment below applies even to the T^O case, we set T = 
for simplicity. For small N, we expand the free energy of 
the ground state up to fourth order in the charged pion 
condensate N = 2\/27r+7r~, 

N 2 N 4 
O = Q + a(n,(i e )— +(3(fi,fi e ) — , 

where Slo is the free energy at N = 0. It can be shown 
that the second order coefficient a has the same struc- 
ture as the 7r + propagator, with Qq replaced by — Mej 
i.e., a(fi, n e ) = jT v + n - (— /i e ; /i, fi e ). Provided that the 
transition to the tt c condensed phase is of second or- 
der, the critical condition is given by the linearized gap 
equation, a(fi, /i e ) = 0, that determines the line of the 
critical electric chemical potential nl(fJ.) in the (n,fx e ) 
plane. From the equality of the pion propagator and 
a noted above, one immediately obtains the threshold 
= ±fj, e for positive (negative) fj, e . This is because 
r^+Tr- (—He, Me) = is also satisfied on the critical line 
//g (//) , and this solution corresponds to the negative (pos- 
itive) root of the BS when fi e > (fj, e < 0). 

The above argument holds either with or without the 
neutrality condition. The neutrality condition determines 
the neutrality line /i" eut (/i) in the (fx, fx e ) plane. At 
to = 10 keV, Me (a*) an d Me eut (M) intersect at the point 
(M, Me) — (324, 10) MeV. This fact is somehow implicit in 
Fig. Q] where the physical quantities only along /Xg eut (//) 
is sketched. We see from the figure that, as soon as /i 
exceeds the value of the vacuum constituent quark mass 
M = 310 MeV, Me becomes finite due to the appearance 
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FIG. 1: The constituent quark mass M, the pion condensate 
N, and meson masses as a function of {i at T = in the neutral 
phase for a toy value of the current quark mass m = 10 keV. 



of quark Fermi surfaces. Accordingly, the mass degener- 
acy in the pion sector gets lifted. 

When the critical point m = 324MeV(= 
Mci) hereafter 2 ) is reached, [i e hits M n - from be- 
low, resulting in the gapless excitation in the 7r~ mode. 
Thereafter, tt~ mode in the N = phase changes 
into the exact Nambu-Goldstone boson associated with 
U T3 (1) breaking; the excitation stays gapless due to 
the mi xing of (a, 7r + ,7r~) modes in the BS equation for 
iV^O [If. 

We also notice that the increase of fi e is somehow 
quenched as soon as the system experiences the tran- 
sition to the pion condensed phase. This is reasonable 
because the tt~ condensate tries to accommodate the d 
quark Fermi surface that helps cancel the positive charge 



2 The critical point [i c i for the transition from the chiral sym- 
metry broken phase to the pion condensed phase corresponds 
to the critical baryon density ng ~ 0.1 n s , with n s being the 
nuclear saturation density 0.16fm -3 . The electron fraction is 
Y e = n e /ng ~ 0.0006. These values are unconventionally small, 
which is simply attributed to the fact that we are working near 
the unphysical chiral limit, m = lOkeV. 
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arising from u quarks; then less electrons are needed. The 
pion condensed phase continues up to /x = 330MeV(= 
/i C 2, hereafter) where the system encounters the first or- 
der phase transition into the almost chirally symmet- 
ric quark matter phase; the phase is characterized by 
M ~ m. Since we are working near the chiral limit 
to = lOkcV, the mass splitting in the sigma and neutral 
pion sector is tiny. In fact M- RQ is only slightly larger than 
M CT ; the difference is less than 1 MeV which is actually 
invisible in the figure. 

A few remarks are in order with regards to the na- 
ture of pion condensed phase obtained here: (i) Having 
a look at the fermionic spectrum, we can see that it is 
of gapless type [23|]. (ii) An investigation of the spectral 
functions in the bosonic sector leads us to conclude that 
it is the BEC of the real du bound state (tt~) rather 
than the BCS-like pionic supcrfluid. Let us first address 
point (ii). It is enough for our present purpose to focus 
on the onset of the pion condensed phase. At the onset 
A* = Met) the 7T~ mode becomes gapless, which means the 
inverse tt~ propagator is zero at to = 0. On the other 
hand, looking at the spectral density, the imaginary part 
of the 7r~ propagator, it is easy to see that the thresh- 
old for the (du) continuum at rest (Q = 0) is located at 
Wth = 2M - /!«(= (M -n- Me/3) + (M + n- 2/ie/a)). 
Then, (u> = 0) in the tt~ propagator is realized as an iso- 
lated bound state pole if 2 M > /x e , while in the opposite 
case it is realized as a soft mode peak whose width goes 
to zero only when its momentum Q goes to zero. Actu- 
ally the former is realized here, so the system is the BEC; 
the isospin chemical potential |/x/| = p, e is too small to 
realize the BCS-like superfluid state [H, [2^. Let us now 
turn to point (i) . The dispersion relations for quasiquarks 
in the pion condensed phase are given by 

Eu(p) = ^(VMl+pl+tle^+N? — (fl— Me/6), 

E 3 (p) = v /( v /M2+p2- Ale /2)2+iV2_ fa _ Me / 6 ), 

( 3 ) 

where the subscript u (d) represents the quasiquark 
which has a large u-quark (d-quark) content at large p. 
Since we have already noted that M > /i e /2 in the BEC, 
each quasiquark energy has a minimum at p = 0. Then 
we see that both quasiquarks have a blocking region in 
p-space, defined by the condition E~,g(p) < 0; the con- 
stituent u and d quarks are, accordingly, accumulated in 
those regions. From (i) and (ii), we may conclude that 
the system is the BEC of the real du bound state (tt~) 
formed in the charge neutral background characterized 
by the Fermi seas of u, d quarks, and electrons. 
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FIG. 2: Phase diagram of neutral matter in (fi, m-n) plane. 



the point of view of calculations but somehow unphysi- 
cal. In this section we investigate the role of the current 
quark mass in 7r c condensation. In this analysis we fix the 
cutoff A and the coupling G to the values specified above 
and we treat m as a free parameter. As a consequence, 
the pion mass at fx = T = 0, m„ in the following, is a free 
parameter as well. In Fig. [2] we report the phase diagram 
in (fi, to^ plane in the neutral case. The solid line rep- 
resents the border between the two regions where chiral 
symmetry is broken and restored. The bold dot is the 
critical endpoint of the first order transition. The shaded 
region indicates the region where 7r c condensation occurs. 
In the chiral limit (m OT = 0) our results are in good agree- 
ment with those obtained in Ref. (23[. As we discussed 
in the previous section, there exist two critical values of 
the quark chemical potential, p, c \ and p C 2, corresponding 
to the onset and vanishing of ir c condensation, respec- 
tively. Increasing the current quark mass results in the 
shrinking of the shaded region till the point p c \ = p C 2 
for m£ ~ 9 MeV, corresponding to a current quark mass 
of to ~ 10 keV. Hence we have shown that the gapless 
7r c condensation is extremely fragile with respect to the 
symmetry breaking effect of the current quark mass; this 
can be understood by observing the fact that the increase 
of quark mass leads to a drastic magnification of the vac- 
uum pion mass at small m because of to w cx \fm, while 
the change in quark mass in keV scale brings about no 
significant modification in \x e . 



VI. THE PHASE DIAGRAM IN THE 
PHYSICAL CASE 



V. THE ROLE OF THE CURRENT QUARK 
MASS 

In the past years, a large number of the calculations 
dealing with the issue of 7r c condensation have been per- 
formed in the chiral limit. This treatment is simple from 



As a final investigation, in Fig. Owe draw the phase di- 
agram of quark matter in the (/x, fj- e ) plane in the physical 
limit where the current quark mass is tuned to to = 5.5 
MeV. At each value of (/x, p e ) we compute the chiral 
and pion condensates by minimization of the thermody- 
namical potential. The solid line represents the first or- 
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FIG. 3: Phase diagram in (fi, fj, e ) plane at m — 5.5 MeV. 

der transition from the ir c condensed phase to the chiral 
symmetry broken phase without the tt c condensate. The 
bold dot is the critical endpoint for the first order tran- 
sition, after which the second order transition sets in. 
The dashed line indicates the first order transition be- 
tween the two regions where chiral symmetry is broken 
and restored, respectively. The dot-dashed line is the neu- 
trality line /i" out = fJ-e(fJ-) which is obtained by requiring 
the global electrical neutrality condition, dCl/dfj, e = 0. 
The neutrality line manifestly shows the impossibility of 
finding a tt c condensate in a physical situation. 

VII. CONCLUSIONS 



up to \l ~ 300 MeV. Furthermore, we have found that 
the possibility of 7r c condensation is ruled out in neutral 
quark matter with a current quark mass larger than the 
order of 10 keV. We have given a natural explanation for 
this fact by using the threshold criterion noted above. 
As a further result we have computed the phase diagram 
in the (//, /x e ) plane at the physical quark mass. We find 
that the transition to the tt c condensed phase is of second 
order up to \i m 300 MeV, after which it turns into a first 
order one. By superimposing the neutrality line onto the 
phase diagram, we have shown that neutral quark matter 
never meets the regions of the nonvanishing ir c conden- 
sate. Thus we conclude that in neutral quark matter at 
the physical quark mass, the NJL model does not allow 
charged pion condensation in the ground state. 

In the present analyses, we have focused on the modifi- 
cation of pion propagation due to the character change of 
the QCD ground state with respect to the finite baryon 
(and isospin) density environment. It should be noted, 
however, that in our restricted treatment, the contribu- 
tion to finite baryon density comes only from constituent 
quarks and not from nuclcons. The existence of nuclc- 
ons, the bound states of constituent quarks, should make 
S-wave 7r c condensation even less favorable. This is be- 
cause the Tomozawa- Weinberg interaction between the 
nucleon and the pion is isospin odd, giving rise to re- 
pulsive pion self-energy at the physical situation where 
the neutron density is larger than the proton density, 
i.e., (n^n) > (p^p) pSO) . It would be interesting, though 
challenging, to investigate pion condensation by taking 
into account nucleon degrees of freedom within the NJL 
model. The prescription made in [3l| may be useful. 



In this article we have studied the phenomenology of 
7r c condensation in two-flavor neutral quark matter us- 
ing a Nambu-Jona-Lasinio model of QCD. In particular, 
we have clarified the role of the quark mass m and the 
electric chemical potential fj, e . Our central results is that 
the threshold to the ir c condensed phase, p e = M x -, 
persists even for a composite pionic gas in the presence 
of dynamical quarks provided the transition to the ir c 
condensed phase is of second order, which is indeed true 
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